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We investigate an extended version of the quantum Ising model which includes beyond-nearest 
neighbour interactions and an additional site-dependent longitudinal magnetic field. Treating the 
interaction exactly and using perturbation theory in the longitudinal field, we calculate the energy 
spectrum and find that the presence of beyond-nearest-neighbour interactions enhances the minimum 
gap between the ground state and the first excited state, irrespective of the nature of decay of these 
interactions along the chain. The longitudinal field adds a correction to this gap that is independent 
of the number of qubits. We discuss the application of our model to implementing specific instances 
of 3-satisfiability problems (Exact Cover) and make a connection to a chain of flux qubits. 
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I. INTRODUCTION 

One of the main motivations for developing scalable 
quantum processors is the realization that carefully con- 
structed quantum algorithms running on such processors 
can solve certain problems that cannot be solved by clas- 
sical computers^. In practice, however, the implemen- 
tation of quantum algorithms using actual qubits - for 
example solid-state qubits such as spin qubits^ or super- 
conducting circuits^ - will be hampered by the presence 
of decoherence, which destroys the interference properties 
on which succesful execution of these algorithms relies. 
In order to try to avoid decoherence effects, Farhi et al. 
proposed in 2001 a method of implementing quantum al- 
gorithms which relies on the adiabatic theorem®. The 
basic idea behind this method, now commonly known as 
adiabatic quantum computation, is to construct a Hamil- 
tonian iffinai whose (unknown) ground state encodes the 
solution to the problem to be solved. By initializing the 
qubit system in the known ground state of a well-chosen 
initial Hamiltonian if initial and letting initial evolve suf- 
ficiently slowly into -ff nn ai, e.g. using 

H(t) = ( 1 — — "\ //initial + j, -^filial, (1) 

the system will end up at t = T in the ground state of 
^finai- Reading out this state then provides the sought - 
for solution of the problem. 

Since the original proposal of Farhi et al., who numeri- 
cally investigated the required running time of adiabatic 
evolution towards a system whose ground state encodes 
the solution of Exact Cover 3 (a NP-completeS' problem 
which belongs to the class of 3-satisfiability problems), 
a lot of research has been done on adiabatic quantum 
computation. The efficiency of adiabatic quantum com- 
putation has been investigated for well-known spin mod- 
els su ch a s the quantum Ising model and the Heisenberg 
model^ and the occupation of the ground state has been 
predicted to be quite robust against decoherence (at suf- 
ficiently low temperatures and for weak coupling of the 



qubit to the environment jSHUl The relation between adi- 
abatic quantum evolution and quantu m phase transitions 
is an ongoing topic of researchl- 11 * 12 !. Also, recently the 
statistics and scaling of energy gaps between the ground 
state and excited states - which form the limiting factor 
for the efficiency of adiabatic quantum computation as 
well as the role played by the choice of -fffinai have been 
investigated^^. 

So far, theoretical proposals for the implementation 
of adiabatic quantum computing have considered m ostly 
generic spin models, such as the quantum Ising mode l 6 ^ 16 l 
These models by themselves cannot be used to encode 
the solution to one of the hard NP-complete problems 
and also in general do not directly correspond to experi- 
mental qubit systems, which are often described by more 
complex versions of these spin models^. 

In this paper we present a first step towards bridging 
the gap between well-understood generic spin models and 
the more complex spin models required for implementing 
adiabatic quantum computing protocols. Specifically, we 
consider an extended version of the quantum Ising model, 
which differs from the standard quantum Ising model in 
two ways: it allows not only for nearest-neighbour, but 
also for next-nearest-neighbour and beyond-next-nearest- 
neighbour interactions, and it includes an additional 
site-dependent longitudinal magnetic field. Building on 
a general exact expression for the energy spectrum of 
this extended quantum Ising model for uniform beyond- 
nearest-neighbour interactions wc include the longitudi- 
nal field using perturbation theory. We analyze the scal- 
ing of the energy gap between the ground state and first 
excited states as a function of the number of beyond- 
nearest-neighbour interactions M and show that the gap 
increases with M both for interactions which decay lin- 
early as a function of the distance between two qubits 
along the chain and interactions which decay exponen- 
tially. We then investigate how our model could be used 
to implement and test particular instances of Exact Cover 
3 that are characterized by limited distance between the 
bits in each clause (corresponding to the maximum num- 
ber of beyond- nearest-neighbours taken into account). 
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We estimate that the probability of errors to occur is 
reasonably low (< 10 %) provided enough neighbours are 
taken into account and enough clauses are defined. Wc 
also discuss the feasibility and prospects of implement- 
ing the extended quantum Ising model using a chain of 
superconducting flux qubits. 

The paper is organized as follows. In Sec. [Il] the prob- 
lem of Exact Cover 3 is introduced, followed by the pre- 
sentation of our model in Sec. IIII1 Sec. QV] contains the 
main calculations: the diagonalization of the quantum 
Ising model with beyond-nearest-neighbour interactions 
(Sec. IV A ) , the modification of the resulting energy spec- 
trum by an additional longitudinal field (Sec. IV B I, and 
the scaling behavior of the gap (Sec. IV C ). We then dis- 
cuss in Sec. |V| how our model can be used to test particu- 
lar instances of Exact Cover 3 and make a connection to 
chains of superconducting flux qubits. Conclusions are 
presented in Sec. VI 



II. 



EXACT COVER 3 



Exact Cover 3 belongs to the class of satisfiability prob- 
lems that are NP-complete^. The problem is the follow- 
ing: a string of N bits X\ . . . x^, which take values or 1, 
has to satisfy M constraints called clauses. Each clause 
applies to three bits, say x a , xp and x^ with a, /3, 7 
€ {1, . . . , N} and is satisfied if and only if one of the bits 
is 1 and the other two are 0: 



X Q/ I X I X -y 1 . 



(2) 



The solution of Exact Cover 3, if it exists, consists of an 
assignment of the bits which satisfies all of the M clauses. 
Of particular interest are instances of Exact Cover with 
a unique solution^. 

In the literature two types of -fffi na i (see Eq. 0) have 
been considered for Exact Cover problems. One involves 
three- qubit interactions^ and the other two-qubit inter- 
actionaHEl] j n both cases, -fffi na i is constructed by asso- 
ciating each violated clause with a fixed energy penalty 
using the "cost function" Y^ all clauscB {x a +x /3 + x 1 - l) 2 . 
In case of two-qubit interactions, i?fi na i is obtained by re- 
placing x a by the Ising variables er" = 1 — 2x a — ±1 and 
substituting cr" by the Pauli operators <r". This yields 
(omitting an irrelevant constant)^ 



N N 
Fflnal = 4 E M <*P<4 ~ ^ E N «°% 



(3) 



a=l 



where er" denotes the Pauli matrix for qubit a (omitting 
the hat), N a represents the number of clauses involving 
qubit a and M a p denotes the number of clauses which 
involve both qubit a and qubit /3. 

From an experimental point of view, both the Hamil- 
tonian involving three-qubit interactions from ReP and 
the Hamiltonian ([3]) are not easy to realize. In existing 
solid-state qubit systems so far three-qubit interactions 



have not been realized yet. Common potentially scal- 
able qubit systems, e.g. electron spin qubit:P or super- 
conducting qubits^' involve two-qubit interactions whose 
strength is a function of the distance between the qubits 
rather than dictated by the clauses (as in Hamiltonian 
([3])). All in all, theoretical predictions of Exact Cover 3 
and other 3-satisfiability problems still seem somewhat 
removed from experimental verification. The aim of this 
paper is to provide a first step towards bridging this gap 
between theory and experiment, by analyzing the model 
Hamiltonian (an extended version of the quantum Ising 
model) that is introduced in the next section. 



III. MODEL 



Our starting point is the time-dependent spin-chain 
Hamiltonian 



E 

i=l 



N IN N 

H(t) = f{t) A < + 9(t) £ ./, .a> 

i=l \i,j=l 

'(4) 

Here N denotes the number of qubits along the chain, 
Jij represents the Ising interaction between qubit i and 
qubit j, A denotes a transverse magnetic field and hi is 
a site-dependent longitudinal field. The functions f(t) 
and g(t) model the time evolution from t = to t = T. 
In most of this paper we choose f(t) = C — t/T, with 
C a constant, and g(t) = t/T. When we deviate from 
this time dependence, this is indicated in the text. For 
any < t < T the instantaneous Hamiltonian Q repre- 
sents an extended quantum Ising model which includes a 
site-dependent longitudinal field hi and whose interaction 
term not only allows for nearest-neighbour interaction 
Ji^i+i but also for next-nearest-neighbour interactions 
and beyond. A similar model was recently considered by 
Amin and Choi^l, who investigated the occurrence of first 
order quantum phase transitions in an inhomogeneous 
version of the Hamiltonian Q. The standard quantum 

Ising model iJi sing = J JV=i + ^Y,^=i a l witn 
fixed nearest-neighbour interactions is a well-known and 
exactly solvable spin modeP^ which has been studied ex- 
tensively for more than 50 years. In the context of adia- 
batic quantum computing, Murg and CiracP have inves- 
tigated adiabatic evolution in the quantum Ising model 
using the ratio A/ J as the time-dependent parameter 
and calculated the excitation probability from the ground 
state to higher-energy states. More recently the robust- 
ness of adiabatic passage against noise was studied^. 

Eq. Q represents a chain of qubits which initially at 
t = is described by the Hamiltonian 



N 



initial 



(5) 
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and has evolved after time t — T into the Hamiltonian form as (omitting an overall minus-sign )p2l 



N 



final 



»=i 



N 



N 



(6) 



i=l 



with f = CA and T = (C - 1)A. The initial Hamilto- 
nian Eq. (|5l describes a chain of spins in a magnetic field 
directed along the z-axis, whose ferromagnetic ground 
state consists of a large superposition of states. The final 
Hamiltonian Eq. ^ reduces to the Hamiltonian ^ for 
r = (i.e. C = 1) and thus encodes the solution of a par- 
ticular instance of Exact Cover 3 if hi (Jy) is interpreted 
as the number of clauses containing bit i (both bit i and 
bit j). For Jij = J and hi = h site-independent, the 
Hamiltonian Eq. ^ reduces to a quantum Ising model in 
an additional uniform longitudinal field h. Using pertur- 
bation theory in h, the ground state of this Hamiltonian 
and the scaling behavior of the gap has been studied in 
RefP. 



tn 



JY 



JY 



»=1 i=l 



Jo 



Jl 



JY 



X]( C 1 ~ C *X C l+2 +C i+2 ) 



^(4_ CiV )(4 + C2 )( e ^ + i) 

^(cj v „ 1 -c Ar _ 1 )( C t 1 +c 1 )(e- L + l). 



(8) 



Here c\ and Cj denote fermionic raising and lowering op- 
erators and L = J2j=i c ) c j as m 

RefP. The last three 
terms are absent in case of periodic boundary conditions, 
and can be neglected for N ^> 1. Diagonalizing ^ using 
Pfeuty's method^ yields: 



Ho 

r 



^AkvUk - J E Afc ' 



(9) 



IV. CALCULATIONS 

In this section we first diagonalize the Hamiltonian ([6| 
in absence of the longitudinal field hi , assuming uniform 
nearest-neighbour interactions and uniform next-nearest- 
neighbour three-qubit interactions^. From the energy 
spectrum we calculate the energy gap between the ground 
state and the first exited state and derive the condition 
for this gap to be minimal. In Sec. |IVB| we then include 
the site-independent longitudinal field hi and calculate 
the corrections to the energy spectrum due to this field 
up to second order in perturbation theory. In Sec. |IV C| 
we use this modified spectrum to analyze the scaling be- 
havior of the gap as a function of the coupling strengths 
A,. 



with the fermionic operators 

N 



E 



4>ki + ^ki\ , ( <\>k% -ipki\ + 
Ci + I I c! 



(10) 



For N even, the sums over k in Eq. ([9| run from —N/2 
to (N - 2)/2. For N odd the sums run from (1 - N)/2 
to (N - l)/2. Defining \j = Jj/(2T), the functions <j) kl 
and ipki are given by 



>ki = 



ipki = - 



N blU V JV 



for k > 
cos(^) for k < 



(11a) 



+ 



uU 1+ lCOS \n~J~ 2COS V*r). 

( _ fc)i )(nb) 



A 



Ai sin 



2vrfc 

1\T 



A 2 sin( — 1 



A. Diagonalization 



and the energy eigenvalues are 



Our starting point is the Hamiltonian: 



JV 



JV 

=1 



(7) 



H originates from the Hamiltonian iffi na i [Eq. ^} by 
taking hi = Mi, defining J\ = J,:,i+i, J 2 = Ji,i+2 Vi, 
taking Jjj = otherwise and adding the third-qubit in- 
teraction in the J 2 -term. By applying a Jordan- Wigner 
transformation to Hq, Eq. ([7| can be rewritten in bilinear 



A? 



/ 47rfc 



1 + Aicos^— J-A 2 cos v N 

\ • f 27Tk \ \ 
Ai sm — — — A 2 sin 



/47rfc\l 2 
\ N J \~N~) 

= l + A? + ^ + 2A 1 (l-A 2 )cos — 



/47rfcA 
-2A 2 cos( — 1 



(12) 



The diagonalization procedure leading to the energy 
spectrum (12 1 can be generalized for higher Aj, j = 



4 



3, 4, . . . (i.e. including interactions between qubits that 
are farther apart) by adding terms 



(-l) j+1 Xj cos 



/ 2irjk\ 



and (-l) i+1 A 



. / 2ir jk 
Sm [-N- 



inside the square brackets in the expressions for tp^ and 
A|. However, one should keep the number M of neigh- 
bours included much smaller than N/2, in order to be 
able to neglect the boundary terms in Eq. (fsh. The full 
expression for the energy spectrum now readipSl 



M 



{2irjk\ 



1 + £(-1)^008 1 — . 

B-^'W 2 ?) 



Note that when iV is even, k = —N/2, and J2jLi A 



(13) 



= 1 



we obtain from Eq. ( 13 ) that A/. = 0. In that case the 



groundstate and first excited state are degenerate and 
the concept of adiabatic transport no longer applied 
After the inclusion of the longitudinal field (see Eq. ([6 
term in the Hamiltonian Hq this degeneracy is lifted, as 
we show below in Section |IVB| For now, however, 
restrict ourselves to the case N odd. From Eq. ([9| 
obtain that the ground state energy of the system is given 
by 



we 
we 



E n 



(14) 



The energy difference between the ground state and the 
excited single fc-fermion state is A& for k running from 
— (N - l)/2 to (N - l)/2. The minimum gap A cg , min 
between the ground state and the first excited state is 
then derived by minimizing A& with respect to a contin- 
uous variable k. Including up to next-nearest-neighbour 
interactions, the minimum is attained at k = ±(JV— l)/2 
for A 2 < 1 and at k = for A2 > 1. In what follows, we 
proceed with the case A2 < 1, which corresponds to the 
most physical situation. Including nearest-neighbours 
only, the energy gap A cg is then given by: 



Aeg=I\ l 



/7T 



2tt 



JV»1 



A? + Al-2Ai(l- A 2 )cos(— ) -2A 2 cos 

(Ai + A 2 - l) 2 + (Ai - AiA 2 + 4A 2 ) (^) 2 . ( 
This gap is minimal for Ai + A 2 = 1, and then equals 
A, 




FIG. 1: (color online) Plot of the single- fermion energy levels 
A fc [Eq. fl2| |] as a function of Ai for N = 51, A 2 = 0.1 and k 
ranging from to ± 25. The lowest graph represents k — ±25. 
From top to bottom (on the right hand side of the plot) the 
graphs are ordered from lowest (k = 0) to highest ranging 
(k — ±25) values. 



The minimum gap A cg . m j n (M ) is thus inversely propor- 
tional to the number of qubits^, irrespective of the num- 
ber of nearest-neighbour interactions M. From Eq. (17) 



we obtain that adding beyond-nearest-neighbour interac- 
tions increases the minimum gap by 



A eg , min (M) - A eg , min (M - 1) = 



In Fig. 1 the energy levels for the single-fermion excited 
states (12) are plotted as a function of Ai for fixed A 2 . 
The minimum gap indeed occurs for Ai ± A 2 = 1. 




B. Perturbation theory in the longitudinal field 

We now apply perturbation theory to calculate the ef- 
fect of the longitudinal field hi [the third term in Eq. ^ 
,which we denote by -f/3] on the level spectrum (13) . In 
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terms of the fermionic operators 77^ and rjk in Eq. (10), 
3 is written as 



H3 

r 



. mm = ^ VAi - AiA 2 ± 4A 2 = ^- (1 ± A 2 ). (16) 



N 



N 



For coupling with more neighbours while keeping N 3> 
M, Eq. ( 16 1 generalizes to 



N h 

En-i 

i=i 



E 

k 



(-D*'(4+«) 
(4 +vk), 



with 



M M 

with E a j = L 

j=i j=i 



hi 



rk = ^(-l) N '< 



-1 

ik 1 



(19) 



(20) 
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Nf the total number of c-fermions and 4> i} ^ given by 
Eq. (26 1. We denote the vacuum state by |0) and the 



state with one 77-fermion by \9 p ). For more rj fermions 
we use more indices, for example \9 Pt q, r ). To first order 
in hi, the correction of the energy of the vacuum state is 
given by 



seW = (o\h 3 \o) 



0. 



(21) 



Analogously, the first-order correction of the energy cor- 
responding to state \9 m ) is given by® 



&E$ = (O m \H 3 \0 m ) 

= -r(e m \J2n (vi + vk) \9 r , 



and are given by (for arbitrary N) 



= (l Sni. Si 

2 2 -■■ 



>ki 



^k = 



1 + Ai cos 
Ai sin 



2?rfc 



/4Trk\ 



2nk 



-A 2 sm^ — 



h i(-k) 

A fc 



hk 
A fe 

(26) 



Finally, we note that Eq. ([20]) can also be used to cal- 
culate r,t if the longitudinal field is site-dependent. For 
example, if hi has a given statistical distribution, Eq. ( 20 ) 



can be used to calculate the corresponding distribution 
and average of r k . 



. k^m 



r k (0m\0k, m ) +r m (e m \0) 



0. 



This line of reasoning can be extended to all odd-order 
corrections of the energy levels, see Appendix |Aj The 
lowest nonzero correction to the energy spectrum is thus 
the second-order correction. For the ground state |0) this 
is given by: 



5E. 



(2) 



E t n (vl + vi) W 



k 



En — Ei. 



^ A fe 



(22) 



To second order the corrections to the energies of the 
single fermion states \0 m ) are given by: 



se^ = r 2 



\(0 m \r mV l\0)\ 2 



Eo 



\{e m \r k y k \0n h k)\ 2 



A r 



fe v K 



(23) 



Higher-order corrections and a discussion of the validity 
of perturbation theory are given in Appendix [A] For a 
site-independent longitudinal field hi = h Vi we obtain 
from Eq. pOl: 



2T 
h_ 

4r 
h_ 

2V 



w tan(^) for k > 

2_ 

N 

2_ 

N 



for k = 
for k < 



(24) 



Here we have used the inverse of the functions 
ipki [Eqns. Ill])], that are defined by 



>ki 



and 



C. Scaling behavior of the gap 



Using Eqns. (22 1 and (23) we now investigate the 



(2) 

second-order correction of the gap SA eg \ m between the 
ground state and the m th single-fermion state: 



SA^ = SE^ 2 



5E t 



(2) 




tan 2 (^) 


rn 


> 


A m 


1 

4A 


rn 


= 




rn 


< 



(27) 



From Sec. |IV A| we know that the minimum gap in the 

absence of the longitudinal field occurs for m = 

and Ej=i = 1< The longitudinal field lifts the degen- 



racy of the ±(^^)-fermion levels. Using Eq. (27) we 



then obtain for the corresponding corrections of the gap 



(at J2f=i Aj = 1): 



SA 



(2) 



h 2 tan (I - ^ t 
FN A n-i 

2 

4h 2 N 1 



SA 



(2) 



Trr 

h 2 



An-i 

2 

1 



FN Ai 



< SA 



(2) 



(28a) 



(28b) 



Substituting Eq. (13) for k = into Eq. (28b), the 

second-order correction to the minimum gap is thus given 
by: 



SA. 



(2) 



(29) 



^<t>ik^ki = X^V'ifcVfct = 1 Vi, 

k k 



(25) Combining Eqns. (17) and ([29]) we then obtain for the 



minimum gap, up to second order in the longitudinal 
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FIG. 2\ (color online) Plot of the minimum gap A eg , m in 
[Eq. ( 30 1] in units of T as a function of M for h/T = 
0.1 and N = 51. The blue dashed line represents linear 
(Aj = const/j) and the red solid line represents exponential 
(Xj = const ■ exp (— j)) decay of the interaction strength with 
distance between the bits. The constant is determined by the 
condition Ylj=i = !• 



field h, 



eg, man 



(M) 



3=1 



(30) 



with ^2j—i Xj = 1. Eq. (30) is the main result of our 
paper. We see that the presence of the longitudinal field 
leads to an increase of the gap by a factor that is in- 
dependent of the number of qubits N. Although the 
/i-dependent correction term Eq. ( 29 ) decreases when 



adding beyond-nearest-neighbour interactions (scaling as 

- ^S min (Af - i) ~ -tfAWfrr) < o) 

the minimum gap itself increases with M since the first, 
unperturbed, term increases with Am- Fig. IV C depicts 



*eg,mxn 



in units of Y as a function of the number of neigh- 
bouring interactions M. We consider both interactions 
which decay linearly as a function of the distance between 
two qubits along the chain and interactions which decay 
exponentially. For both types of decay (although more 
strongly for linear than exponential decay) the minimum 
gap is enhanced by including coupling with more neigh- 
bours. Since the required running time of an adiabatic 
algorithm is inversely proportional to the square of the 
energy garP, the enhancement of the gap implies that 
a qubit system with beyond-nearest-neighbour interac- 
tions may be advantageous for implementing adiabatic 
quantum algorithms. 



V. RELATION TO EXACT COVER 3 

In this section we investigate how the extended quan- 
tum Ising model from the previous section with uni- 



form beyond-nearest-neighbour interactions can be used 
to simulate specific instances of Exact Cover 3. In par- 
ticular, we numerically estimate the probability for ob- 
taining the correct solution of instances of Exact Cover 3 
that are characterized by a maximum distance between 
bits in a clause as a function of the number M of beyond- 
nearest-neighbour interactions that are included. At the 
end of the section we make a connection to an actual 
experimental system (a chain of flux qubits). 

We consider a system of N bits with coupling to M 
nearest neighbours. Since each bit is coupled to M bits 
on either side, the maximum distance between two bits 
that appear in the same clause is 1M (see also Eq. (|3|). 
Out of the set of all possible clauses, we restrict ourselves 
to clauses that satisfy this property, i.e. contain bits 
which are at most a distance 2M along the chain apart. 
We refer to this subset as the set of "restricted clauses". 
Clauses in which no maximum distance between the bits 
is defined, are called "unrestricted clauses". For clauses 
{a, /?, 7} with a < f3 < 7 we only allow restricted clauses 
that satisfy 7 - a < 2M or (a + N) - < 2M or 
(/3 + N) — 7 < 2M and in order to avoid boundary ef- 
fects we consider a cyclic chain of bits. The question 
that we raise is the following: what is the probability 
that a particular solution of Exact Cover 3, represented 
by a random bit chain (in which each bit independently 
has a probability p for taking value 1) can be recon- 
structed using K restricted clauses? In our simulations 
the restricted clauses for given K are selected from a 
large group of unrestricted clauses that are uniformly dis- 
tributed over the N bits in the qubit chairPn For some 
limits the answer to this question is obvious: in partic- 
ular, for p close to or 1 there are too many bits with 
the same value, and no restricted set of clauses, or in- 
deed unrestricted set of clauses, can be found that gives 
a uniquely satisfying assignment. 

We now describe our simulations. For specific values 
of N, M, p, and K we generate 100 random bit strings^. 
For each of these we randomly generate 500 clauses (or 
more if needed) that are satisfied by the given bit string. 
Out of these we randomly select K clauses that satisfy 
the M-nearest-neighbour restriction condition. Then we 
check whether the original bit string is the unique solu- 
tion of the restricted set of clauses. If each bit is covered 
by at least one clause, the solution is usually unique. 
Out of 100 runs we deduce a probability pe that an er- 
ror will occur, meaning that at least one other bit string 
also satisfies the assignment of K clauses. Our goal is 
to determine at which point pe makes a transition from 
large (pe ~ 1) to small (pe ~ 0) as a function of the 
parameters p and M. In the table below, we keep N 
and K fixed at N — 12 and K = 20, and calculate pe 
as a function of K for various probabilities p. The ratio 
K/N is chosen large enough that once the bits allow a 
uniquely satisfying assignment of restricted clauses for a 
given p, this is also formed in the simulations with high 
probability. 

From Tabic [j] we see that for all values of p a transition 
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M 


1 


2 


3 


4 


5 


6 


7 


8 


Pe for p = 0.2 


1.00 


0.64 


0.36 


0.22 


0.14 


0.14 


0.10 


0.09 


Pe for p = 0.3 


1.00 


0.43 


0.22 


0.11 


0.10 


0.10 


0.09 


0.05 


Pe for p = 0.4 


1.00 


0.51 


0.26 


0.18 


0.16 


0.10 


0.10 


0.09 


Pe for p = 0.5 


1.00 


0.67 


0.34 


0.23 


0.22 


0.18 


0.17 


0.13 



TABLE I: The probability pe for errors to occur as a function 
of the number M of beyond-nearest-neighbour interactions for 
p = 0.2,0.3,0.4, and 0.5. 



1.0 
0.8 
0.6 
0.4 
0.2 



ft, 




whose final ground state at time T encodes the solution 
of the Exact Cover problem. In their general form, the 
Hamiltonians Q and Q are still far from experimen- 
tal realization. The closest analogy between Q and an 
actual experimental qubit system is probably a chain of 
coupled flux qubits, which can be described by the Hamil- 
toniarPMl 



N 



N 



-Hfluxchain = ^ Ai<7* +^ Ki^ +1 a\a 



i+1 



N 
i=l 



(31) 



Here e« denotes the magnetic energy, Aj the tunnel 
coupling energy of individual qubits and K i i+ i ~ 
Mi^ + ilplp +1 with M iyi+ i the mutual inductance between 
the persistent currents 7* and of qubits i and i + 
respectively. In practice, local (in situ) tuning of the pa- 
rameters A,, Ki t i + i and ej, as required by adiabatic evo- 
lution of Eq. Q, is challenging, but promising progress 



is being mad 



This gives hope for achieving con- 



trol over individual qubits in chain and array-like qubit 
geometries and thereby brings experimental simulation 
of adiabatic quantum algorithms such as Exact Cover 3 
closer. 



VI. CONCLUSION 



FIG. 3: (color online) The probability pe from Table [I] as a 
function of M for p = 0.2 (red, solid), p = 0.3 (blue, dot- 
dashed), p = 0.4 (green, dashed), and p — 0.5 (purple, dot- 
ted). 



from large to small error probability takes places for M 
ranging from 1 to 8. For p = 0.3 the probability of errors 
to occur for a given value of M is smallest, which can be 
explained from the optimal ratio of bitvalues and 1 for 
finding clauses. Fig. |V| shows a graphic representation of 
Table[l] We see that the probability for an error to occur 
during simulation of restricted instances of Exact Cover 
3 decreases approximately exponentially with increasing 
M. 

Our simulations did not take into account the fact that 
the interaction strength A„ between qubits i and i + n - 
which translates into the number of clauses that contain 
both bit i and bit i+n, see Eq. ^ - in practice decreases 
as a function of the distance n between the qubits. In 
order to make a connection to the Exact Cover Hamil- 
tonian ^ we in principle thus need to further restrict 
the sets of "restricted clauses" to sets in which the num- 
ber of clauses containing nearby qubits is larger than the 
number of clauses containing bits that are farther apart. 
Although we did not investigate this in depth, back-of- 
the-envelope estimates indicate that in order to reach the 
same success probability a factor of 2 more clauses need 
to be included in the simulations. 

In order to simulate Exact Cover 3 in practice, one 
needs to use adiabatic evolution of the Hamiltonian (ffl , 



In conclusion, we have calculated the energy level spec- 
trum of the quantum Ising model in the presence of uni- 
form beyond-nearest-neighbour interactions and an addi- 
tional longitudinal field. We found that the gap between 
the ground state and the lowest-lying excited state in- 
creases with increasing number M of neighbouring in- 
teractions (provided M remains much smaller than the 
total number of qubits along the chain) and is approxi- 
mately linearly proportional to M for linearly decreasing 
interaction strength between the qubits along the chain. 
The increase of this gap with M, which persists in the 
presence of the additional weak longitudinal field, gives 
hope that the extended quantum Ising model is suitable 
for numerical - and in the future hopefully also experi- 
mental - simulation of quantum algorithms such as Exact 
Cover 3. 
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Appendix A 

In this Appendix we investigate the validity of the per- 
turbative approach that we used in Sec. |IVB] and derive 
a criterion for the application of perturbation theory in 
terms of the number of qubits (N) and the number of 



nearest neighbours included (M). We first demonstrate 
that all odd-order corrections to the single- fermion en- 
ergy levels are zero and then investigate even-order cor- 
rections. Our starting point is the third-order correction, 
which in general form is given by: 



SE^ 



EE- 

q9 Lp r =£p 



9 p \H 3 \9 q )(9 q \H 3 \9 r )(9 r \H 3 \9 p ) 
(E p — E q ) (Ep — E r ) 



(Al) 



It can be seen immediately that the second term is zero 
since it contains the matrix element (9 P \H 3 \9 P ) , which is 



zero (see Eq. (22)). The first term in Eq. (Al) contains 



three separate matrix elements. Let us assume that in the 
first product 9 P is an even-fermion state, which is then 
coupled to odd- fermion states 9 q by H 3 . This in turn 
restricts 9 r to the even fermion subspace. Applying the 
same reasoning to the last inner product we see that for 
this product to be non-zero 9 p should be an odd-fermion 
state. This is a contradiction with our starting assump- 
tion. Hence there is no combination of states that gives 
a non-zero outcome. This observation can be extended 
to all odd-order energy corrections and we conclude that 
these are therefore all zercP^. 

We now investigate even-order corrections to the en- 
ergy of the ground state and single-fermion states and use 
these to derive a criterion for the validity of perturbation 
theory. Our starting point is the second-order correc- 
tion of the ground state, which can be found by inserting 
Eqns. (pi) and p4b into Eq. p2l. This yields: 



6E, 



(2) 



h 2 



2NT 




(A2) 

After defining M = Y]j—i jAj, rewriting the sum using 
m = N — 2k and expanding around m = (which yields 
the largest contribution to the sum), Eq. ( A2) reduces to 



8E, 



, 2 N-l , /2JVV 
(2) _ h ST" 1 + \™) 

2NT ^ M(W-) 

m— 1 \ N / 

m odd 

h 2 N 2 
v -0.068- 



TM 



(A3) 



Analogously we obtain for the fourth-order correction: 

= r£Q;) 2 (!) 2 ^(A l + A„r' 



k.l 

4 h 4 N 5 



N-l N-l _A _9 

m n 



E E 



7T 5 (TM) 3 ^ ^ m + n 

v / m— 1 n— 1 
m odd n odd 

w 0.0071 ^^. (A4) 
(TM) 3 

We are interested in the energy gap between the ground 
state and the first excited states (N — 2k <C N), which 
in the absence of the longitudinal field h is given by 



- E^ 



4 0) = TA fc = ^TM(N - 2k), (A5) 



and reduces to Eq. (17) for the lowest- lying excited state 
(given by k = (N — l)/2). The second-order corrections 
to this gap between the ground state and the lowest-lying 
excited states (N — 2k <C N) are given by (see Eq. (27)) 



SE, 



(2) 



SE, 



(2) 



2T 



{r k f 



4 h 2 



N 



" A fe tt3 r M(N-2k)' 



(A6) 



For the general (2n) -order correction to the gap we find 



?i2ni\r3rt— 1 

SE^-SE^k^S . 

k (rM) 2 "- 1 



(A7) 



It follows from Eq. ( A7 ) that perturbation theory is valid 
for 



h M 

f < AT3/2' 



Fig. IV C in the main text shows a plot of A 



eg,min : 



(A8) 



which 



is directly proportional to M, as a function of the number 
of neighbouring interactions M. For exponential decay 
the value of M converges for M large (but still much 
smaller than TV), but for linear coupling M scales almost 
linear with M. For a given number of qubits N and linear 
decay of interaction strength along the chain adding more 
nearest-neighbours thus enhances the range of validity of 
perturbation theory. 
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